Extended matter coupled to BF theory 



Winston J. Fairbairn * 
Laboratoire de Physique - CNRS-UMR 5672 
ENS Lyon 
46, Allee d'ltalie 
F-69364 Lyon, EU 

Alejandro Perez ' 
Centre de Physique Theorique x - CNRS-UMR 6207 
Luminy - Case 907 
F-1S288 Marseille, EU 

o 

Q 1 (Dated: February 2, 2008) 

Recently, a topological field theory of membrane-matter coupled to BF theory in arbitrary space- 
O I time dimensions was proposed In this paper, we discuss various aspects of the four-dimensional 

D ■ theory. Firstly, we study classical solutions leading to an interpretation of the theory in terms of 

strings propagating on a flat spacetime. We also show that the general classical solutions of the 
theory are in one-to-one correspondence with solutions of Einstein's equations in the presence of 
distributional matter (cosmic strings). Secondly, we quantize the theory and present, in particu- 
lar, a prescription to regularize the physical inner product of the canonical theory. We show how 
the resulting transition amplitudes are dual to evaluations of Feynman diagrams coupled to three- 
dimensional quantum gravity. Finally, we remove the regulator by proving the topological invariance 
of the transition amplitudes. 
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J> ' Based on the seminal results [|| of 2 + 1 gravity coupled to point sources, recent developments [||, [f| 

*0 ! in the non-perturbativc approach to 2 + 1 quantum gravity have led to a clear understanding of quantum 

field theory on a three-dimensional quantum geometrical background spacetime. The idea is to first couple 
free point particles to the gravitational field before going through the second quantization process. In this 
approach, particles become local conical defects of spacetime curvature and their momenta are recasted as 
holonomies of the gravitational connection around their worldlines. It follows that momenta become group 
valued leading to an effective notion of non-commutative spacetime coordinates. The Feynman diagrams of 
| such theories are related via a duality transformation to spinfoam models. 

All though conceptually very deep, these results remain three-dimensional. The next step is to probe all 
possible extensions of these ideas to higher dimensions. Two ideas have recently been put forward. The first 
is to consider that fundamental matter is indeed pointlikc and study the coupling of worldlines to gravity by 
using the Cartan geometric framework Q of the McDowell-Mansouri formulation of gravity as a de-Sitter 
gauge theory Q . The second is to generalize the description of matter as topological defects of spacetime 
curvature to higher dimensions. This naturally leads to matter excitations supported by co-dimension two 
membranes @. Before studying the coupling of such sources to quantum gravity, one can consider, 
as a first step, the BF theory framework as an immediate generalization of the topological character of 
three-dimensional gravity to higher dimensions. 

This paper is dedicated to the second approach, namely the coupling of string-like sources to BF theory 
in four dimensions. The starting point is the action written in [l[ generating a theory of flat connections 
except at the location of two-dimensional surfaces, where the curvature picks up a singularity, or in other 
words, where the gauge degrees of freedom become dynamical. The goal of the paper is a two-fold. Firstly, 
acquire a physical intuition of the algebraic fields involved in the theory which generalize the position 
and momentum Poincare coordinates of the particle in three-dimensions. Secondly, provide a complete 
background independent quantization of the theory in four dimensions, following the work done in 
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The organization of the paper is as follows. In section II, we study some classical solutions guided by the 
three-dimensional example. We show that some specific solutions lead to the interpretation of rigid strings 
propagating on a flat spacctime. More generally, we prove that the solutions of the theory arc in one-to-one 
correspondence with distributional solutions of general relativity. In section III, we propose a prescription 
for computing the physical inner product of the theory. This leads us to an interesting duality between 
the obtained transition amplitudes and Feynman diagrams coupled to three-dimensional gravity. We finally 
prove in section IV that the transition amplitudes only depend on the topology of the canonical manifold 
and of the spin network graphs. 

II. CLASSICAL THEORY 

A. Action principle and classical symmetries 

Let G be a Lie group with Lie algebra q equipped with an A<i(G)-invariant, non degenerate bilinear form 
noted 'tr' (e.g. the Killing form if G is semi-simple). Consider the principal bundle V with G as structure 
group and as base manifold ad + l dimensional, compact, connected, oriented differential manifold M. We 
will assume that V is trivial, all though it is not essential, and chose once and for all a global trivialising 
section. We will be interested in the following first order action principle, describing the interaction between 
closed membrane-like sources and BF theory [l[ : 

S[A, B; q,p] = S BF [A, B] - f tr{B + d A q)p). (1) 

The action of free BF theory in d + 1 dimensions is given by 

S B f[A,B] = - [ tr{BAF[A}). (2) 

K J M 

Here, B is a g-valued (d — l)-form on M, F is the curvature of a g-valued one-form A, which is the pull-back 
to M by the global trivializing section of a connection on V, and k <E K is a coupling constant. 

In the coupling term, W is the (d — l)-branc worldshcet defined by the embedding (j> : E C — > M, 
d A is the covariant derivative with respect to the connection A, q is a g-valued (d — 2)-form on W and p 
is a g-valued function on W . The physical meaning of the matter variables p and q will be discussed in 
the following section. Essentially, p is the momentum density of the brane and q is the first integral of the 
(d — l)-volume element; the integral of a line and surface element in in three and four dimensions {d = 2, 3) 
respectively. 

The equations of motion governing the dynamics of the theory are those of a topological field theory: 



F[A] = up S w (3) 

d A B = n\p, q]6yp (4) 

<P*{B + d A q) = (5) 

d A p\w = 0- (6) 



Here, Syp- is a distributional two-form, also called current, which has support on the worldsheet W . It is 
defined such that for all (d — l)-form a, JL, a = J M (a A 8-p-). The symbol (f>* denotes the pull-back of forms 
on W by the embedding map (f>. 

We can readily see that the above action describes a theory of local conical defects along brane-likc (d—1)- 
submanifolds of M through the first equation. The second states that the obstruction to the vanishing of the 
torsion is measured by the commutator of p and q. The third equation is crucial. It relates the background 
field B to the dynamics of the brane. For instance, this equation describes the motion of a particle's position 
in 3d gravity [2j. The last states that the momentum density is covariantly conserved along the worlsheet. 
It is in fact a simple consequence of equation ([3]) together with the Bianchi identity d A F = 0. We will see 
how this is a sign of the reducibility of the constraints generated by the theory. 

The total action is invariant under the following (pull back to M of) vertical automorphisms of V, 

V.g e C°°{M, G), B ^ B = gBg" 1 (7) 
A i-> A = gAg^ 1 + gdg' 1 
P !-> 9P9 1 

q !-> gqg' 1 



3 



and the 'topological', or reducible transformations 

Vr7ef2 d - 2 (M,£0, 

where tt p (M, g) is the space of g- valued p-forms on M. 



B i ► B + d A ri (8) 

p i— > p 
9 !-* 9 - ?7 



B. Physical interpretation: the flat solution 

In this section, we discuss some particular solutions of the theory leading to an interpretation of matter 
propagating on flat backgrounds. We discuss the d = 2 and d = 3 cases where the gauge degrees of 
freedom of BF theory become dynamical along one dimensional worldlincs and two-dimensional worldshects 
respectively. 

1. The point particle in 2 + 1 dimensions 

We now restrict our attention to the d = 2 case with structure group the isometry group G = SO (77) 
of the diagonal form r/ of a three-dimensional metric on M; 77 = (a 2 , +, +) with a = in respectively 

Riemannian (G = SO(3)) and Lorcntzian (G = SO(l, 2)) signatures. We denote (ir, V v ) the vector (adjoint) 
representation of 50(77) = R{ J a }o=o,i,2, i-e., V v = M 3 and V v = K 1,2 in Riemannian and Lorentzian signatures 
respectively. The bilinear form 'tr' is defined such that tr( J a Jb) = jVab- In this case, the free BF action @ 
describes the dynamics of three-dimensional general relativity, where the B field plays the role of the triad e. 
The matter excitations are 0-branes, that is, particles and the worldshcct W reduces to a one-dimensional 
worldlinc that we will note 7. The degrees of freedom of the particle are encoded in the algebraic variables 
q and p which are both so (revalued functions with support on the world-line 7. 

Firstly, we consider the open subset U of M constructed as follows. Consider the three-ball B 3 centered 
on a point xq of the worldline 7 and call x and y the two punctures dB 3 n 7. Pick two non intersecting 
paths 71 and 72 on dB 3 both connecting x to y. The open region bounded by the portion of dB 3 contained 
between the two paths and the two arbitrary non intersecting disks contained in B 3 and bounded by the 
loops 771 and 772 defines the open subset U C M. 

Next, we define the coordinate function X : M — > V v mapping spacctime into the 'internal space' 50(77) 
isomorphic, as a vector space, to its vector representation space Vq. The coordinates are chosen to be 
centered around a point x in M traversed by the worldlinc; X(x) = 0. Associated to the coordinate 
function X, there is a natural solution to the equations of motion d3j), (|4j), (J6]), (JSj) in C/ 

e = dX = 5 (9) 

A = 
q=-X\ 1 
p = constant, 

where 8 is the unit of End(T p M, V v ), 5(v) = v forall v in T p M and all p in U. The field configuration e = 5 
(together with the A = solution) provides a natural notion of flat Riemannian or Minkowskian spacetime 
geometry via its relation to the spacetime metric g = 2tr(e <E> e). This flat background is defined in terms 
of a special gauge (notice that one can make e equal to zero by transformation of the form ©)■ From now 
on, we will call such gauge a flat gauge. The solution for q is obtained through the equation §5§ relating 
the background geometry to the geometry of the worldline. Here, we can readely see that q represents the 
particle's position X, first integral of the line element defined by the background geometry e. Below we 
show that equation (J3J) forces the worldline to be a straight line. Finally, p = constant trivially satisfies 
the conservation equation ([B]). In fact, the curvature equation of motion ([3]) constrains p to remain in a 
fixed adjoint orbit so we can introduce a constant m 6 R* + such that p = mv with v S 20(77) such that 
tre 2 = —a 2 . Consequently, p satisfies the mass shell constraints p 2 := tr(p 2 ) = —a 2 m 2 and acquires the 
interpretation of the particle's momentum. 

We can now relate the position q and momentum p. independent in the first order formulation, by virtue 
of d3]). Indeed, the chosen flat geometry solution e = S, A = leads to a everywhere vanishing torsion d^e. 
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Hence, the commutator [p, q] = X x p, where x denotes the usual cross product on V v , vanishes on the 
worldlinc. This vanishing of the relativistic angular momentum (which is conserved by virtue of equation 
(J5J)) implies, together with the flatness of the background fields, that the the worldline 7 of the particle 
defines a straight line passing through the origin and tangent to its momentum p. Equivalently, we can 
think of the momentum p as Hodge dual to a bivector *p, in which case the worldlinc is normal to the plane 
defined by *p. 

Note that translating 7 off the origin, which requires the introduction of spacetime torsion, can be achieved 
by the gauge transformation q — > q + C with C = constant which leaves all the other fields invariant. In 
this way we conclude that the previous solution of our theory can be (locally) interpreted as the particle 
following a geodesic of flat spacetime. 

More formally, we can also recover the action of a test particle in flat spacetime by simply 'switching 
off' the interaction of the particle whith gravity. This can be achieved by evaluating the action |T]) on the 
flat solution and neglecting the interactions between geometry and matter, namely the equations of motion 
linking the background fields to the matter degrees of freedom (e.g. e 7^ dX). This formal manipulation 
leads to the following Hamilton function 

S[p, X, N] = f tx(pX) + N(p 2 - m 2 ), (10) 
which is the standart first order action for a relativistic spinless particle. 



2. The string in 3 + 1 dimensions 

We now focus on the four dimensional (d = 3) extension of the above considerations. Here again we 
consider the isometry group G = SO(t]) of a given four dimensional metric structure 77 = (a 2 , +, +, +), in 
which case the value a = 1 leads to the Riemannian group G = SO (4), while a = i encodes a Lorentzian 
signature G = SO(l,3). As in three dimensions, we denote (n, V v ), with V v = E{ej}j, I = 0, ...,3, the 
vector representation of so(n) = R{ J a b}a,b=o,...,3- Finally, we choose the bilinear form 'tr' such that, forall 
a, b in so(n), it is associated to the trace tr(a&) = harjb IJ in the vector representation. We are using the 
notation a IJ = a ab n( J a b) IJ '■= ot ab J J a ( for the matrix elements of the image of an element a £ so(n) in 
End(y i; ) under the vector representation. The dynamics of the theory is governed by the action ((T|) where 
the matter excitations are string-like and the worldshect W is now a two-dimensional submanifold of the 
four dimensional space time manifold M. The string degrees of freedom are described by an so (revalued 
one-form q and an so(?7)-valued function p living on the world-sheet W . 

As before, we construct an open subset U C M by cutting out a section of the four-ball B A , and define 
the coordinate function X : M — > V V: centered around a point x in M R W . Consider the following field 
configurations which define a flat solution to the equations of motion ([3]), ©, © in U .: 

B = *(eAe), with e = dX = 8 (11) 
A = 

q = -*XdX 
p = constant, 

where the star V is the Hodge operator * : il p (V v ) ~ > il 4 ~ p (V v ) acting on the internal space; (*a)u = 
ie 7 J <L ctKL, with the totally antisymmetric tensor e normalized such that e 0123 = +1. 

The solutions B = *(<5A<5) (A = 0), leads to a natural notion of flat Riemannian or Minkowski background 
geometry through the standard construction of a metric out of B when B is a simple bivector; B = *(e A e) 
with e ~ S. We can readily see that the q one-form is the first integral of the area element defined by the 
background field B. As in 3d, the equations of motion constrain p to remain in a fixed adjoint orbit so that 
we can introduce a constant r £ R* + such that p = tv and v £ 50(77) has a fixed norm; trv 2 = —a 2 . We 
call r the string tension, or mass per unit length, and p the momentum density which satisfies a generalized 
mass shell constraint. 

This momentum density p is related to the q field by analysis of equation The solution (B = *5/\5, A = 
0) has zero torsion d^-B- Accordingly, the commutator [p, q] = [*XdX,p] vanishes on the worldsheet. This 
leads to the constraint X I pu = 0. Putting everything together, we see that the flat solution in the open 
subset U leads to the picture of a locally flat worldshect (a locally straight, rigid string) in flat spacetime, 
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dual 1 , as a two-surface, to the momentum density bivector p (if p is simple, namely if it defines a two-plane) . 
If we consider more general solutions admitting torsion, the plane can be translated off the origin. Indeed, 
the equation (|5|) determines the field q in terms of the geometry of the B field up to the addition of an exact 
one-form j3 = da which encodes the translational information. For instance, the translation X X + C of 
the plane yields q i— > q — *CdX and consequently corresponds to a function a defined by da = — * CdX. 
This potential is in turn determined by the torsion T = dAB of the B field via the equation More 
general solutions can be found for arbitrary a's, as discussed below. 

Following the same path as in the case of the particle case, we can 'turn off' the interaction between the 
topological BF background and the string by evaluating the action on the flat solution (this implies, here 
again, that we ignore the equations of motion of the coupled theory, i.e. the relation between the matter 
and geometrical degrees of freedom). We obtain the following Hamilton function 



up to a constant. This is the Polyakov action on a non trivial background with metric G M „ = and 
antisymmetric field b = v G 50(77). 

Now, the previous action leads to trivial equations of motion that are satisfied by arbitrary X (because 
p = constant and so the Lagrangian is a total differential). This is to be expected, from the string theory 
viewpoint, this would be a charged string moving on a constant potential, so the field strength is zero. This 
seems in sharp contrast to the particle case where the effective action leads to straight line solutions. Here 
any string motion is allowed; however, from the point of view of the full theory, all these possibilities are 
pure gauge. The reason for this is that in 2 + 1 dimensions the flat gauge condition e = 6 fixes the freedom 
© up to a global translation, and hence gauge considerations are not necessary in interpreting the effective 
action. In the string case, B = *(5 A 6) partially fixes the gauge; the remaining freedom being encoded in 
T) = da for any a. 



The above discussion shows that particular solutions of the theory in a particular open subset lead to 
the standard propagation of matter degrees of freedom on a flat (or degenerate) background spacctimc. In 
fact, we can go further in the physical interpretation by considering other solutions, defined everywhere, 
which are in one-to-one correspondence with solutions of four dimensional general relativity in the presence 
of distributional matter. These solutions are called cosmic strings. 



It is well known that the metric associated to a massive and spinning particle coupled to three-dimensional 
gravity is that of a locally flat spinning cone. The lift of this solution to 3 + 1 dimensions corresponds to 
a spacetime around an infinitely thin and long straight string (see for instance ff5j and references therein). 
Let us endow our spacetime manifold M with a Ricmannian structure (M, g) and let x G M label a 
point traversed by the string. We can choose as a basis of the tangent space T X M the coordinate basis 
{dt, d r , d z } associated to local cylindrical coordinates such that the string is lying along the z axis 
and goes through the origin. The embedding of the string is given by (f>(t,z) = (t, 0,0, 2). Let t and s 
respectively denote the mass and intrinsic (spacetime) spin per unit length of the string. Note that r is the 
string tension. Solving Einstein's field equations for a such stationary string carrying the above mass and 
spin distribution produces a two-parameter (r, s) family of solutions described by the following line element 
written in the specified cylindrical coordinates 




(12) 



C. Geometrical interpretation: cosmic strings and topological defects 



1. Cosmic strings 



ds 2 = g^dx^ 1 ® dx v 

= a 2 (dt + 0dtp) 2 + dr 2 + (I - a) 2 r 2 dip 2 + dz 2 , 



(13) 



Note that this is exactly the same result than the one obtained for the point particle, if we think of the 3d momentum as 
Hodge dual to a bivector. 
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where (3 = AGs and a = (1 — 4GV). G is the Newton constant. In fact this family of metrics is the general 
solution to Einstein's equations describing a spacetime outside any matter distribution in a bounded region 
of the plane (r, ip) and having a cylindrical symmetry. Exploiting the absence of structure along the z axis, 
by simply suppressing the z direction, reduces the theory to that of a point particle coupled to gravity in 
2 + 1 dimensions, where the location of the particle is given by the point where the string punctures the 
z = plane. We will come across a such duality again in the quantization process of the next sections. The 
dual co-frame for the above metric is written 

e° = dt + l3dip (14) 
e = cos ipdr — ar sin ipd(p 
e 2 = sin ipdr + ar cos (pd(p 
e 3 = dz, 

such that ds 2 = e 1 ®e J 'rjjj.li we assume that the connection A associated to the above metric is Ricmannian, 
it is straight-forward to calculate its components by exploiting Cartan's first structure equation (d^e = 0). 
The result reads 

A = Aj/ aijdx 11 = AGt a 12 dip, (15) 

where {er/j}/,,/ is a basis of Q 2 (V V ) ~ so(r/). 

Using the distributional identity ddip = 2irS 2 (r)dxdy (x = rcos(p, y = r sirup, and dxdy is a wedge 
product), it is immediate to compute the torsion T = T°eo and curvature F = F 12 oyi of the cosmic string 
induced metric : 

T° = 8nGs S 2 (r) dxdy, F 12 = SttGt 8 2 {r) dxdy. (16) 

These equations state that the torsion and curvature associated to the cosmic string solution are zero 
everywhere except when the radial coordinate r vanishes, i.e. at the location of the string worlshcet lying 
in the z — t plane. If we now focus on the spinless cosmic string case s = 0, we can establish a one-to-one 
correspondence between the above solutions of general relativity and the following solutions of BF theory 
coupled to string sources: 

B al = sin ipdrdz + ar cos tpdtpdz, B 02 = — (cos tpdzdr — ar sin ipdzdp) 

B 03 = ardrdip, B 12 = -a 2 dzdt 

B 13 = — a 2 (sin ipdtdr + ar cos tpdtdip), B 23 = a 2 (cos ipdtdr ~ ar sin tpdtdip), 



A 12 = 4Grd(p, (17) 
q 12 = a 2 (zdt-tdz), p 12 = t, 

where only the non vanishing components have been written and the coupling constant n in (TTJ) has been 
set to 8irG. 

In this way, solutions of our theory are in one-to-one correspondence to solutions of Einstein's equations. 
The converse is obviously not true as our model does not allow for physical local excitations such us 
gravitational waves. However, augmenting the action |T]) with a Plebanski term constraining the B field to 
be simple, would lead to the full Einstein equations in the presence of distributional matter, 

euK L e J A F KL = 8ttGt e IJKL e J J{\ L 8 W , (18) 
where J(^ L = 5^ 5 2 \ starting from the theory considered in this paper. 



2. Many-strings-solution 

One can also construct a many string solution by 'superimposing' solutions of the previous kind at 
different locations. Here we explicitly show this for two strings. We do this as the example will illustrate 
the geometric meaning of torsion in our model. Assume that we have two worlsheets W\ and #2 respectively 
traversing the points pi and p 2 - We will work with two open patches E/j C M , i = 1,2, such that p\ and 
p 2 both belong to the overlap U\ (~l U 2 - The cylindrical coordinates (ti,n,<pi, z{) associated to the charts 
(Ui C M, X? : Ui — > IR 4 ) are chosen such that the strings lie along the z axis, are separated by a distance 
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Xq in the x-direction, and are such that ri(j>i) = 0. The coordinate transform occurring in the overlap 
U\ n U 2 is immediate; it yields U = t, x 2 = x\ + xq yi = y and Zi = z, for i = 1, 2. The two embeddings are 
given consequently by 4>i (t, z) = (t, 0, 0, z) and fait, z) = (t, x , 0, z). Our notations are such that a field <j> 
expressed in the coordinate system associated to the open subset Ui is noted <j)u i . 

Our strategy to construct the two-string-solution is the following. We need to realize the fact that, 
regarded from a particular coordinate frame, one of the two strings is translated off the origin. We will 
choose to observe the translation of W 2 from the coordinate frame 1. Now, the study of the flat solution 
discussed in the previous section has showed that translations of the worlsheet are related to the torsion T of 
the B field. In particular, we know how to recognize a translation of the form X — > X + C, with C = x§e\. 
It corresponds to a torsion of the form T = K[p,da], with da — — * CdX. Hence, the two-string-solution 
is based on the tetrad field which leads to the desired value of the B field torsion taking into account the 
separation of the two worldsheets. For simplicity, here we assume that the two strings are parallel, hence 
that they have same momentum density 

Pm = PU 2 = rcr 12, (19) 

and accordingly create the same curvature singularity in both coordinate frames 1 and 2. The associated 
connection yields 

Au^AGrdtpiaw, Vi = l,2. (20) 
The dual co-frame ejj i — e\j. ® eju i is defined by the following components 

e%. = dt 

e\j. = cos (fidri — ari sin ifidifii 
efj = sin (fidri + (an cos tfi + S i2 -^-Tx )dipi 

47T 

efj. — dz. 

By integrating the B = *e A e solution with e given by (|2"Tj) , we can now calculate the q field, up to the 
addition of an exact form f3 — da 

qu i = a 2 (zdt — tdz) oyi + da\ J ajj. (22) 

The potential a is derived from the equation of motion relating the commutator of p and q to the B 
field torsion three-form T = dAB = *dAe A e + *e A dA& ■ 

Tu z = S i2 - k tx S(r)(dx 2 dy 2 dz cr i + o -2 dt dx 2 dy 2 013). (23) 

This torsion indeed corresponds to a two-string-solution since it yields the desired value — * CdX for the 
form da, 

da t = 5 l2 i x (dz 002 + o 2 dt <r 23 ). (24) 

One can add more than one string in a similar fashion, leading to multiple cosmic string solutions. It is 
interesting to notice that torsion of the mutiplc string solution is related to the distance xq separating the 
world sheets. Of course this is a distance defined in the flat-gauge where B = *6 A 5. This concludes our 
discussion on the physical aspects of the action (JTJ) of string-like sources coupled to BF theory. We now 
turn toward the quantization of the theory. 



(21) 



III. QUANTUM THEORY 



For the entire quantization process to be well defined, we will restrict our attention to the case where the 
symmetry group G is compact. For instance, we can think of G as being SO (4). We will also concentrate on 
the four-dimensional theory and set the coupling constant k to one. Also, to rely on the canonical analysis 
performed in [l| , we will work with a slightly different theory where the momentum p is replaced by the 
string field A E C°°(W ,G). This new field enters the action only through the conjugation rAd\(v) of a 
fixed unit element v in g, and the theory is consequently defined by the action (JJ) with p set to rAuA -1 . 
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The field A transforms as A — > gX under gauge transformations of the type ||7J) and the theory acquires a 
new invariance under the subgroup H C G generated by v. The link between the two theories is established 
by the fact that, as remarked before, the equation of motion F = p8yp- implies that p remains in the same 
conjugacy class along the worldsheet. Here, we choose to label the class by tv and to consider A as dynamical 
field instead of p. 



As a preliminary step, we assume that the spacetime manifold M is diffeomorphic to the canonical split 
R x E, where R represents time and E is the canonical spatial hypersurface. The intersection of E with the 
string worldsheet W forms a one dimensional manifold that we will assume to be closed 2 . We choose 
local coordinates (t, x a ) for which E is given as the hypersurface {£ = 0}. By definition, x a , a = 1, 2, 3, are 
local coordinates on E. We also choose local coordinates (t, s) on the 2-dimensional world-sheet W , where 
s £ [0, 2tt] is a coordinate along the one-dimensional string & '. We will note x :9 . = <j> |s the embedding of the 

string 5? in E. We pick a basis {X,} i=1 dim{a) 01 the- real Lie algebra g, raise and lower indices with the 

inner product 'tr', and define structure constants by [.Xfj-Xj] = Next, we choose a polarization on 

the phase space such that the degrees of freedom are encoded in the configuration variable (A, A) £ A x A 
defined by the couples formed by (the pull-back to E of) connections and string momenta. 

The canonical analysis of the coupled action |T]) shows that the legendre transform from configuration 
space to phase space is singular : the system is constrained. Essentially 3 , the constraints are first class and 
are given by the following set of equations : 



Here, Ef = e abc Bn, c is the momentum canonically conjugate to Af, irf = d s x%,pi is conjugate to q and 
satisfies D a irf = 0, where pt = tr(JQp) denote the components of the Lie algebra element p in the chosen 
basis of g. The symbols D and F denote respectively the covariant derivative and curvature of the spatial 
connection A. 

The first constraint (|25[). the Gauss law, generates kinematical gauge transformations while the second 
(1261) . the curvature constraint, contains the dynamical data of the theory. To quantize the theory, one 
can follow Dirac's program of quantization of constrained systems which consists in first quantizing the 
system before imposing the constraints at the quantum level. The idea is to construct an algebra 21 of basic 
obscrvables, that is, simple phase space functions which admit unambiguous quantum analogues, which is 
then represented unitarely, as an involutive and unital *-algebra of abstract operators, on an unphysical or 
auxiliary Hilbcrt space H. Since the classical constraints are simple functionals of the basic observables, 
they can be unambiguously quantized, that is, promoted to self-adjoint operators on 7i. The kernel of these 
constraint operators are spanned by the physical states of the theory. 

The structure of the constraint algebra enables us to solve the constraints in different steps. One can 
first solve the Gauss law to obtain a quantum kinematical setting. Then, impose the curvature constraint 
on the kinematical states to fully solve the dynamical sector of the theory. In [l[ , the kinematical subset of 
the constraints is solved and a kinematical Hilbert space 7i kin solution to the quantum Gauss law is defined. 
We first review the kinematical setting of [lj before exploring the dynamical sector of the theory. 



The Hilbert space W kin of solutions to the Gauss constraint is spanned by so-called string spin network 
states. String spin network states are the gauge invariant elements of the auxiliary Hilbert space Ti of 
cylindrical functions which is constructed as follows. 



A. Canonical setting 



d := njr; ■ « o 

H? := e abc F lbc ~^5y w 0. 



(25) 
(26) 



B. Quantum kinematics : the Gauss law 



2 In fact if S is compact the equations of motion <£3j implies that the string must be closed (or have zero tension). 

3 See the original work [lj for a detailed canonical analysis. 
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1. Auxiliary Hilbert space TL 

Firstly, we define the canonical BF states. Let r C E denote an open graph, that is, a collection of 
one dimensional oriented sub-manifolds 4 of E called edges er, meeting if at all only on their endpoints 
called vertices vp. The vertices forming the boundary of a given edge er are called the source s(er) and 
target t(er) vertices depending on the orientation of the edge. We will call n = rip the cardinality of 
the set of edges {er} of T. Let <j> : G Xn — > C denote a continuous complex valued function on G xn 
and A(er) = g er — Pexp(J A) denote the holonomy of the connection A along the edge er £ T. The 
cylindrical function associated to the graph T and to the function <f> is a complex valued map ^r^ : A — > C 
defined by : 

*r, [A] = ^(e r ),...,^(e?)), (27) 

forall A in A. The space of such functions is an abelian *-algebra denoted CylsF.r, where the ^-structure 
is simply given by complex conjugation on C. The algebra of all cylindrical functions will be called GylsF 
= Ur Cyl_B_F,r- 

Next, wc define string states. Since the configuration variable is a zero-form, we expect to consider wave 
functions associated to points ieS. Accordingly, wc define the *-algebra Cyls of cylindrical functions on 
the space A of A fields as follows. An element &x.f of Cyls is a continuous map &x,f '■ A — > C, where 
X = {xi,...,x n } is a finite set of points in and / : G xn — > C is a complex valued function on the 
Cartesian product G n , defined by 

<S>xj[X}=.f(X(x 1 ),...,X(xn))- (28) 

Both algebras Cyls^ and Cyls, regarded as vector spaces, can be given a pre-Hilbert space structure. Fixing 
a graph r C E with n edges and a set of m points X C E, wc define the scalar products respectively on 
CyWr and Cyls, P as 

<3- r ^,$r,v>> = / M (29) 

JG X " 

and 

<*'xj,*x, B >= f 7.9, (30) 

J G xm 

where the integration over the group is realized through the Haar measure on G. These scalar products can 
be extended to the whole of Cyl^ (resp. Cyls), i.e. to cylindrical functions defined on different graphs 
(resp. set of points), by redefining a larger graph (resp. set of points) containing the two different ones. The 
resulting measure, precisely constructed via projective techniques, is the AL measure. The string Hilbert 
space was in fact introduced by Thiemann as a model for the coupling of Higgs fields to loop quantum 
gravity (l3j via point holonomies. Completing these two pre-Hilbert spaces in the respective norms induced 
by the AL measures, one obtains the BF and string auxiliary Hilbert spaces respectively denoted Hbf and 
TLs- Tensoring the two Hilbert spaces yields the auxiliary Hilbert space TL = Hbf ® TCs of the coupled 
system. 

Using the harmonic analysis on G, one can define an orthonormal basis in H.bf and Tis the elements 
of which are respectively denoted (open) spin networks and n-points spin states. Using the isomorphism 
of Hilbert spaces L 2 (G xn ) ~ er L 2 (G er ), any cylindrical function ^r,<f> in TLbf decomposes according to 
the Peter- Wcyl theorem into the basis of matrix elements of the unitary, irreducible representations of G : 

*r,*[A]= J2 ^Pi....,PnPi[^(ef)]®»-®Pn[^(er)]. (31) 
pi,...,p„ 

where p : G — ► Aut(V p ) denotes the unitary, irreducible representation of G acting on the vector space V p 
and the mode </> Pl; ... jPre := ®" =1 </> Pi is an element of ( Y Pi <S> Y Pi *)®*=i . The functions appearing in the above 
sum are called open spin network states. 



4 More precisely, one usually endows the canonical hypersurface S with a real, analytic structure and restricts the edges to 
be piecewise analytic or semi-analytic manifolds, as a mean to control the intersection points. 
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Equivalcntly, the string cylindrical functions decompose as : 



Pl[X(xi)} ® ... (g) p m [Hx m )}, 



(32) 



pl, — ,pm 

and a given element in the sum is called an n-point spin state. 

2. String spin network states 



One can now compute a unitary action of the gauge group G°°(£, G) on TL by using the transformation 
properties of the holonomies and of the string fields A — > gX under the gauge group and derive the subset of 
G-invariant states, that is, the states solution to the Gauss constraint. A vectorial basis of the vector space 
of g auge invariant states can be constructed, in analogy with 3d quantum gravity coupled to point particles 
[12j, by tensoring the open spin network basis with the n point spin states elements. Such an tensorial 
element is required the following consistency conditions to be G-invariant. 

The graph T of the open spin network has a set of vertices Vr including the points {x\, ...,x n } forming 
the set X. The vertices of T are coloured with a chosen element l v of an orthonormal basis of the vector 
space of intertwining operators 



Home 



e r |t(e r )=t>r 



e r \s(e r )=v T 



(33) 



if the vertex vr is not on the string. If a vertex vr is on the string, it coincides with some point Xk G X . In 



this case, we chose an element l Vt in an orthonormal basis of 



Home 



®v Pfc 



e r \t(e r )=v r 



l ep \s(er)—vr 



(34) 



where V Pfc is the representation space associated to the point Xk ■ 

By finally implementing the invariancc under the sub-group H C G generated by v of the n-point spin 
states by choosing the modes to be H- invariant, one obtains a vectorial basis in the kinematical Hilbert 
space H kin where the inner product is that of (BF and string) cylindrical functions. The elements of this 
basis are called string spin networks states and are of the form (see fig. [1} : 



$r,x[A, A] := (* r ® $x)[A A] = 



g) p er [A(e r )] 0p»[A(x)] 
_e r er xex 



(35) 



where the dot '.' denotes tensor index contraction. 

This concludes the quantum kinematical framework of strings coupled to BF theory performed in [l[ . We 
now solve the curvature constraint and compute the full physical Hilbert space W phyB . 



C. Quantum dynamics: the curvature constraint 



In this section, we explore the dynamics of the theory by constructing the physical Hilbert space 7i phyB 
solution to the last constraint of the system, that is, the curvature or Hamiltonian constraint PU)) . Note 
that the physical states that we construct below are also solutions to the constraints of four dimensional 
quantum gravity coupled to distributional matter, as in the classical case. 

We first underline a crucial property of the curvature constraint of d + 1-dimcnsional BF theory with 
d > 2, namely its reducible character which has to be taken into account during the quantization process. 
We then proceed (as in pjjl) a la Rovelli and Reisenberger [13, HH by building and regularizing a generalized 
projection operator mapping the kinematical states into the kernel of the curvature constraint operator. This 
procedure automatically provides the vector space of solutions with a physical inner product and a Hilbert 
space structure, and leads to an interesting duality with the coupling of Feynman loops to 3d gravity [l7| , 
[14j from the covariant perspective. 
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FIG. 1: A typical string spin network (the string is represented by the bold line). 

1. The reducibility of the curvature constraint 

A naive imposition of the curvature constraint on the kincmatical states leads to severe divergences. This 
is due to the fact that there is a redundancy in the implementation of the constraint; the components of 
the curvature constraint of (d+ l)-dimensional BF theory are not linearly independent, they are said to be 
reducible, if d > 2. The same is true for the theory coupled to sources under study here. As an illustration 
of this fact, let us simply count the degrees of freedom of source free (t = 0) BF theory in d+ 1 dimensions. 

The configuration variable of the theory A l a is a g-valued connection one-form, thus containing d x dim(g) 
independent components for each space point of E. In turn, the number of constraints is given by the dim(g) 
components of the Gauss law (f2"5| plus the d x dim(g) components of the Hamiltonian constraint (f2"6"| for 
each space point x € T,. Hence, we have Nq = (d + 1) x dim(g) constraints per space point. This leads to 
a negative number of degrees of freedom. What is happening 5 ? The point is that the Nc constraints are 
not independent: the Bianchi identity (D^F = d^F + [A, F] = 0, where the superscript ( p ) indicates the 
degree of the form acted upon) imply the reducibility equation 

D a H? = 0. (36) 

In the case where sources are present, the reducibility equation remains valid because the curvature con- 
straint F = pS y together with the Bianchi identity automatically implements the momentum density 
conservation Dp = 0. We will come back to this reducibility of the matter sector of the theory. The system 
is said to be (d — 2)-th stage reducible in the first class curvature constraints. This designation is due to 
the fact that the operator d^ is himself reducible since d^d^ = 0. In turn, d 1 -- 3 * 1 is reducible and so on. 
The chain stops after precisely d — 2 steps since the action of the d^ de-Rham differential operator on 
d-forms is trivial. Accordingly, the Nr = dim(rj) reducibility equations (|36[) imply a linear relation between 
the components of the curvature constraint. The number Ni of independent constraints is thus given by 
Nc — Nr = d x dim(g). Using Nj to count the number of degrees of freedom leads to the correct answer, 
namely zero degrees of freedom for topological BF theory. 

The standard procedure to quantize systems with such reducible constraints consists in selecting a subset 
H \i rr of constraints which are linearly independent and impose solely this subset of constraints on the 
auxiliary states of Ji. 

Keeping this issue in mind, we now proceed to the definition and regularization of the generalized projector 
on the physical states and construct the Hilbert space 7i phya of solutions to all of the constraints of the theory. 



5 We thank Merced Montesinos for pointing out this property of BF theory's constraints. 
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2. Physical projector : formal definition - the particle/ string duality 

We start by introducing the rigging map 

Vphys : Cyl - Cyl* (37) 
* i-> S(H \ irr )^, 

where Cyl* is the (algebraic) dual vector space of Cyl = CylsF ® Cylg. The range of the rigging map ?7 phy3 
formally lies in the kernel of the Hamiltonian constraint of the coupled model. The power of the rigging 
map technology is that it automatically provides the vector space ?7 phy3 (Cyl) = Cyl* h s C Cyl* of solutions 
to the Hamiltonian constraint with a pre-Hilbert space structure encoded in the physical inner product 

< »w(*i).'w(*2) > phys = :=< *i, S(H \„ r )^> 2 >, (38) 

for any two string spin network states \&2 G W kin . The scalar product used in the last equality is the 
kinematical inner product (|29[) . ([BIT]) . The physical Hilbert space 7i phy3 is then obtained by the associated 
Cauchy completion of the quotient of H kin by the Gel'fand ideal defined by the set of zero norm states. 

Accordingly, the construction of the physical inner product can explicitly be achieved if we can rigorously 
make sense of the formal expression 5{H \i rr )- This task is greatly simplified by virtue of the following 
duality. Indeed, we can re express the above formal quantity as follows 

S(H \ irr ) = J] S{H \ irr (x)) = ex P (* J tr ( N A H)j ■ (39) 

Here, Af 3 N is the space of regular g-valued one-forms on £ and T>fi[N] denotes a formal functional measure 
on Af imposing constraints on the test one- form N to remove the redundant delta functions on H. Simply 
plugging in the explicit expression of the exponent in (|3"9"f leads to 

H[N] = I tr(N A H) = f tr(N A F) + f ti(Np) 

= S% d F+part [N,A], (40) 

which, in the case where G = SO(rA, where 77 is a three-dimensional metric, is the action of 3d gravity 
coupled to a (spinless) point particle 01 : the role of the triad is played by N, the mass and the worldlinc 
of the particle are respectively given by the string tension r (hidden in the string variable p = rAd\(v)) 
and '. Finally, the role of the Cartan subalgebra generator J is played by v (also hidden in p). This 
relation is reminiscent to the link between cosmic strings in 4d and point particles in three-dimensional 
gravity discussed in the first sections. More generally, we have in fact the following duality : 

P(tyBF+(d-2)-branes = ^BF+(d-3)-branes^ (41) 

where f2 denotes the (d+ 1-dimensional) no-spin-network vacuum state, Z is the path integral of BF theory 
in d spacetime dimensions and we have introduced the linear form 6 P on Cyl C 21 defined by 

V* € Cyl, P(¥) = < ?/phy3 (f!), ?/phy3 (*) > phys (42) 

= < O ,6(6 \ irr ) * > • 

Furthermore, when d = 3, the formal functional measure T>p[N] introduced above to take into account the 
reducible character of the four-dimensional theory corresponds to the Fadeev-Popov determinant gauge- 
fixing the translational topological symmetry of the 3d theory; the reducibility of the 4d theory is mapped 
via this duality onto the gauge redundancies of the three-dimensional theory. 



More precisely, the linear form P, once normalized by the evaluation P(l), is a state 

P/P(l) : Cyl C 21^ C, 

whose associated GNS construction leads equivalently to the physical Hilbert space 7i p h ys . The associated Gel'fand ideal X 
is immense by virtue of the topological nature of the theory under consideration. Indeed, one can show that any element of 
Cyl based on a contractible graph is equivalent to a complex number. The associated physical representation 7r p h yB : Cyl — ♦ 
End(W p h y s) is defined such that forall cylindrical function ar S Cyl defined on a contractible graph T, 7r phy s(ar [^4]) , I' = 
a 7 [0]*, forall * in H p h ya . 



13 



Now because of the above duality (|41[) . regularizing the formal expression (|39p is, roughly speaking, 
equivalent to regularizing the path integral for 3d gravity coupled to point particles, up to the insertion of 
spin network observables. The physical inner product in our theory will therefore be related to amplitudes 
computed in (27[, although they would have here a quite different physical interpretation. Following [H, 
[l9| , we will regularize the Hamiltonian constraint at the classical level by defining a lattice-like discretization 
of £ and by constructing holonomies around the elementary plaquettes of the discretization as a first order 
approximation of the curvature. 

However, there are two major obstacles to the direct and naive implementation of such a program. The 
first is the reducible character of the curvature constraint and the second is the presence of spin network 
edges ending on the string. We will use the above duality to treat the first issue while the second will be 
dealt with by introducing an appropriate rcgularization scheme. 

3. Physical projector : rcgularization 

Throughout this section, we will concentrate on the definition of the linear form (|^|) evaluated on the 
most general string spin network state \& G 7i kin , since it contains all the necessary information to compute 
transition amplitudes between any two arbitrary elements of the kinematical Hilbert space. We will consider 
a string spin network basis elements ^> of W kin defined on the (open) graph T. The set of end points of the 
graph living on the string £f will be denoted X . 

We follow the natural generalization of the regularization defined in [l2j for 2 + 1 gravity coupled to point 
particles. In order to deal with the curvature singularity at the string location, we thicken the smooth curve 
.y to a torus topology, smooth, non-intersecting tube T v of constant radius r\ > centered on the string 
,9* . The radius 77 is defined in terms of the local arbitrary coordinate system. If the string is disconnected, 
we blow up each string component in a similar fashion. 

Next, we remove the tube T v from the spatial manifold £. We are left with a three-manifold with torus 
boundary £ \ T n noted £ n . For instance, if £ has the topology of S 3 , we know by Heegard's splitting, 
that the resulting manifold has the topology of a solid torus whose boundary surface is the Heegard surface 
defined by the string tube. In this way we construct a new three-manifold with boundary where each 
boundary component is in one-to-one correspondence with a string component and has the topology of a 
torus. Finally, the open graph T is embedded in the bulk manifold and its endpoints lie on the boundary 
torus. 

The next step is to choose a simplicial decomposition 7 of T< v or more generally any cellular decomposition, 
i.e., a homcomorphism <j) : £ I; — > A from our spatial bulk manifold £^ to a cellular complex A. The 
discretized manifold A = A e depends on a parameter e G M. + controlling the characteristic (coordinate) 
'length scale' of the cellular complex. We will see that, by virtue of the three-dimensional equivalence 
between smooth, topological and piecewise-linear (PL) categories, together with the background independent 
nature of our theory, no physical quantities will depend on this extra parameter. We will note A& the fc-cells 
of A. To make contact with the literature, we will in fact work with the dual cellular decomposition A*. 
The dual cellular complex A* is obtained from A by placing a vertex v in the center of each three-cell A3, 
linking adjacent vertices with edges e topologically dual to the two-cells A2 of A, and defining the dual faces 
/, punctured by the one-cells Ai, as closed sequences of dual edges e. The intersection between A* and the 
boundary tube T v induces a closed, oriented (trivalent if A is simplicial) graph which is the one-skeleton of 
the cellular complex dA* = (v, e, /) dual to the cellular decomposition <9A of the 2d boundary T v induced 
by the bulk complex A. We will note T the set of faces / of the cellular pair (A*, dA*) and require that 
each dual face of T admits an orientation (induced by the orientation of £,,) and a distinguished vertex. 

Finally, among all possible cellular decompositions, we select a subscctor of two-complexes which are 
adapted to the graph T. Namely, we consider dual cellular complexes (A*, dA*) whose one-skeletons admit 
the graph I as a subcomplex. In particular, the open edges of T end on the vertices v of the boundary 
two-complex dA*. 

The meaning of the curvature constraint F = p b~y is that the physical states have support on the space 
of connections which are flat everywhere except at the location of the string where they are singular. In 
other words, the holonomy g 7 = ^(7) of an infinitesimal loop 7 circling an empty, simply connected region 



7 Note that in dimension d < 3, each topological cZ-manifold admits a piccewise-lincar-structure (this is the so-called 'triangu- 
lation conjecture'). 
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yields the identity, while the holonomy g 1 circling the string around a point x S 5? is equal to expp(x), the 
image of the fixed group clement u = e TV under the inner automorphism Ad\ : G — ► G; u i— > A(x) u A _1 (x), 
with the string field A evaluated at the point x. The integration over the string field A appearing in the 
computation of the physical inner product then forces the holonomy of the connection around the string to 
lie in the same conjugacy class Cl(u) than the group element u. 

To impose the F = part of the curvature constraint, we will require that the holonomy 

A:T ^ G (43) 
df " 9f = LI A ( e )> 

eCdf 

around all the oriented boundaries of the faces / of T be equal to one 8 . Each such flat connection defines a 
monodromy representation of the fundamental group 7Ti(E^) in G. Concretely, the holonomics arc computed 
by taking the edges in the boundary df of the face / in cyclic order, following the chosen orientation, starting 
from the distinguished vertex. Reversing the orientation maps the associated group element to its inverse. 

It is here crucial to take into account the reducible character of the curvature constraint to avoid di- 
vergences due to redundancies in the implementation of the constraints (i.e. divergences coming from the 
incorrect product of redundant delta functions). As discussed above, the reducibility equation induced 
by the Bianchi identity implies that the components of the curvature are not independent. In the dis- 
cretized framework, we know [l8| that forall set of faces / forming a closed surface >S with the topology of 
a two- sphere, 

II 9f = a ' ( 44 ) 

modulo orientation and some possible conjugations depending on the base points of the holonomies. Ac- 
cordingly, there is, for each three-cell of the dual cellular complex A*, one group element c//, among the 
finite number of group variables attached to the faces bounding the bubble, which is completely determined 
by the others. It follows that imposing gj = 1 on all faces of the cellular complex A* is redundant and 
would create divergences in the computation of the physical inner product. The proper way [22| . [l4j to 
address the reducibility issue, or over determination of the holonomy variables, is to pick a maximal tree 
T of the cellular decomposition A and impose F = only on the faces of A* that are not dual to any one 
simplex contained in T. A tree T of a cellular decomposition A is a sub-complex of the one-skeleton of A 
which never closes to form a loop. A tree T of A is said to be maximal if it is connected and goes through 
all vertices of A. The fact that T is a maximal tree implies that one is only removing redundant flatness 
constraints taking consistently into account the reducibility of the flatness constraints. 

Finally, we need to impose the F = p part of the curvature constraint. The idea is to require that the 
holonomy g 7 around any loop 7 in 9A* based at a point x, belonging to the homology class of loops of 
the boundary torus T v normal to the string 5? (these loops are the ones wrapping around the cycle of the 
torus circling the string, i.e., the non-contractible loops in E^), be equal to the image of the group element 
u under the adjoint automorphism A(x)mA(x) _1 , i.e., belong to Cl(u). Intuitively, this could be achieved 
by picking a finite set {7^}; of such homologous paths all along the tube and imposing g 7i = Ad^iu), 
with the field Aj evaluated at the base point of the holonomy g li . However, here again, care must be taken 
in addressing the reducibility issue induced by the equation D a Hf = 0. In the presence of matter, the 
reducibility implies that the curvature constraint F = p8y together with the Bianchi identity DF = 
induce the momentum density conservation Dp = 0. In our setting, this is reflected in the fact that the 
holonomics g 7l and <7 72 associated to two distinct homologous loops 71 an 72 circling the string satisfy the 
property CZ(<7 71 ) = Cl{g l2 ) on shell. This is due to the Bianchi identity in the interior of the cylindrically 
shaped section of the torus T v bounded by 71 and 72, and the flatness constraint F = imposing the 
holonomies around all the dual faces on the boundary of the cylindrical section to be trivial (see e.g. (|44p). 
Accordingly, imposing g~ ll G Cl{u) naturally implies that g l2 belongs to the conjugacy class labeled by 
u. In other words, choosing one arbitrary closed path circling the string, say 71 based at a point xi, and 
imposing F = p only along that path naturally propagates via the Bianchi identity and the flatness constraint 
and forces the holonomy g 72 around any other homologous loop 72 based at a point X2 to be of the form 



Note that the blow up of the string, reflected here in the presence of a flatness constraint on the boundary torus, gives us 
the opportunity to impose that the connection is flat also on the string. 
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g 12 = huh" 1 6 Cl(u), for some /igG. This shows that imposing F = p more than once, e.g. also around 
72, would lead to divergences which can be traced back to the reducibility of the constraints. 

However, for the prescription to be complete 9 , it is not sufficient to have g l2 in Cl(u); we need to recover 
the fact that the holonomy along the loop 72 is the conjugation of the group element u under the dynamical 
field A evaluated at the base point X2 of the holonomy, namely A2 = \{x2). This suggests an identification 
of the group element h conjugating u with the value of the string field A2, which leads to a relation between 
the holonomy gp along a path (3 connecting the points X\ and xi and the value of the string field A at x\ 
and X2- 

9p = K(p)\(p)~ j 

stating that A is covariantly constant along the string. We have seen that the Bianchi identity together with 
the full curvature constraint induces the momentum density conservation. Our treatment of the reducibility 
issue consists in truncating the curvature constraint, i.e., in imposing F = p only once, and using the Bianchi 
identity supplemented with the momentum conservation Dp = to recover the truncated components of 
the curvature constraint without any loss of information. 

Accordingly, the full prescription is defined via a choice of a closed, oriented path a and a finite set 
C of open, oriented paths j3 in dA*. The closed path a circles the string (it is non-contractible in the 
three-manifold S^). This loop is based at a point x £ X lying on a dual vertex v 6 dA* supporting a spin 
network endpoint. The open paths (3 6 C are defined as follows. Let 7 € dA* be an oriented loop based 
at x, non-homologous to a (along the cycle of T v contractible in E^) and connecting all the spin network 
endpoints Xk G X. Define the open path 7 by erasing the segment of 7 supported by the edge e which is 
such that x = t(e). The paths (3 of C are Id sub-manifolds of 7 each connecting x to a vertex v traversed 
by 7. If the graph T is closed, one reiterates the same prescription simply dropping the requirements on 
the spin network endpoints x^ e X, in particular, the base point x is chosen arbitrarily. We then impose 
g a = expp with p evaluated at the point x, and gp = \ s (p)^t(p)~ l , where x = s{(3), on each open path (3 of 
C. 

To summarize, we choose a regulator H(~ t \ = (T v , (A e , <9A e ), T, a, C) consisting in a thickening T v of the 
string, a cellular decomposition (A e , <9A e ) of the manifold (£ l; , T v ) adapted to the graph T, a maximal tree 
T of A, a closed path a in dA*, and a collection C of open paths (3 in dA* . The associated regularized 
physical scalar product is then given by 

P[*]:= lim P[K in , £ y,y}, (45) 

with 



P[^, ie) ;*]=< 



II S ^f) Y[S(g a cxpp) Yl S(gpX t (p)K(p) 

f£T a p£C 



* >, (46) 



where the product over a is to take into account the possible multiple connected components of the string. 

It is important to point out that, in addition to the expression of the generalized projection above, we can 
use the regularization to give an explicit expression of the regularized constraint corresponding to H [N] in 
equation (|40[) . With the notation introduced so far the regulated quantum curvature constraint becomes 

H V . X [N] = ^[N(x f )g f ] + ^[N(x p )g a cxpp], (47) 

/GA* a 



To understand these last points, consider two loops 71 and 72 belonging to the same homology class circling the string 
at two neighboring points, say x\ and X2- These two loops define a section of the torus T v homcomorphic to a cylinder. 
Suppose that the dual cellular complex A* is such that the cylinder is discretized by a single face with two opposite edges 
glued along an dual edge (3 in dA* connecting x\ and X2- The flatness constraint on the boundary of the cylinder implies 
the following presentation of the cylinder's fundamental polygon: 

9719/3972 1 9/3 1 = 1 ' 

which relates the holonomies g 71 and g 12 by virtue of the Bianchi identity in the interior of the tube. Hence, imposing 
F = p only along one of the two loops, say 71, naturally leads to the constraint <? 72 S Cl(u). Finally, plugging the relation 
9/3 = -^l-^ 1 m * ne va l ue °f the holonomy g 12 leads to the required constraint g 12 = ^uX^ 1 ■ 



16 



where Xf is an arbitrary point in the interior of the face / and x p and arbitrary point on the string dual 
to the loop a (the sum over a is over all the string components). It is easy to check that the regulated 
quantum curvature constraints satisfy off-shell anomaly freeness condition. For instance 

U[g]H v , e [N}tf[g] = H v%e \gNg~\ (48) 

where U[g] is the unitary generator of G-gauge transformations. Therefore the regulator does not break 
the algebraic structure of the classical constraints. The quantization is consistent. The quantum constraint 
operator is defined as the limit where rj and e are taking to zero. Instead of doing this in detail we shall 
simply concentrate on the regulator independence of the physical inner product in the following section. 

The inner product is computed using the AL measure, that is, by Haar integration along all edges of the 
graph (A*,<9A*) 3 T and along all endpoints Xk S X C dA* of the graph T. 

We can now promote the classical delta function on the lattice phase space to a multiplication operator 
on 7i kin by using its expansion in irreducible unitary representations : 

V<?eG, S(g)=Y^dim(p) Xp (g), (49) 
p 

where \p = tr p : G — > C is the character of the representation p. Each Xp is then promoted to a self-adjoint 
Wilson loop operator x p on 7i kin creating loops in the p representation around each plaquette defined by our 
regularization, which is charged for the face bounded by the loop a. To summarize, we have, for each face 
/ of the regularization, a sum over the unitary, irreducible representations pf of G, a weight given by the 
dimension dim(p/) of the representation pf summed over and a loop around the oriented boundary of the 
face in the representation pf. See for instance [T^l, [I{J,[2l} for details. This concludes our regularization 
of the transition amplitudes of string-like sources coupled to four dimensional BF theory. 

Now, the physical inner product that we have constructed above depends manifestly on the regulat- 
ing structure 7Z( v ^y To complete the procedure, we have to calculate the limit in which the regulating 
parameters r\ and e go to zero. 

IV. REGULATOR INDEPENDENCE 

Throughout this section we will suppose that the cellular complex (A,<9A) is simplicial, i.e. a triangu- 

lation of (E^,! 1 ^). We will also make the simplifying assumption that G = SU(2), the unitary irreducible 

j 

representations of which will be noted (n, Vj), with the spin j in N/2. The generalization to arbitrary cellu- 
lar decompositions and arbitrary compact Lie groups can be achieved by using the same techniques that we 
develop below. If A is a now a regular triangulation, it cannot be adpated to any graph T. It can only be 
so for graphs with thrcc-or-four-valent vertices. Now, we can always decompose a n-valent intertwiner, with 
n > 3, into a three- valent intertwiner by using repeatedly the complete reducibility of the tensor product of 
two representations 10 . We will therefore decompose all n-valent interwiners l v with n > 3 into three- valent 
vertices. 

We now show how to remove the regulator TZ^^) from the regularized scalar product (|4"5)) . Instead of 
computing the 77, e — > limit, we demonstrate that the transition amplitudes are in fact independent of the 
regulator. To prove such a statement, we show that the expression (|45[) does not depend on any component 
of the regulator. The transition amplitudes are proven to be invariant under any finite combination of 
elementary moves called regulator moves 31 : ^-(^e) 1— » IZ'^ e ,\, where each regulator move is a combination 
of elementary moves acting on the components of the regulator : 

• Bulk and boundary (adapted) Pachner moves (A,<9A) 1— > (A',<9A'), 

• Elementary maximal tree moves T 1— > T 1 ', 

• Elementary curve moves 7 1— ► 7'. 



One can show that the amplitudes obtained using the three-valent decomposition (where the virtual edges are assumed to 
be real) are identical to the ones obtained using a suitable non-simplicial cellular decomposition adapted to spin networks 
with arbitrary valence vertices. 
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We will see how the invariance under the above moves also implies an invariance under dilata- 
tion/contraction of the string thickening radius r\. 

To conclude on the topological invariance of the amplitudes from the above elementary moves, we will 
furthermore prove that the transition amplitudes are invariant under elementary moves acting on the string 
spin network graph S : T 1— > T' which map ambient isotopic PL-graphs into ambient isotopic PL-graphs. 

We now detail the regulator and graph topological moves. 

A. Elementary regulator moves 

The regulator moves are finite combinations of the following elementary moves acting on the simplical 
complex (A, dA), the maximal tree T and the paths a and /? G C. 

1. Adapted Pachner moves 

The first invariance property that we will need is that under moves acting on the simplicial-pair (A, dA), 
leaving the one-complex T invariant and mapping a T-adapted triangulation into a PL-homeomorphic T- 
adapted simplicial structure. We call these moves adapted Pachner moves. There are two types of moves 
to be considered : the bistellar moves [28|, acting on the bulk triangulation A and leaving the boundary 
simpicial structure unchanged, and the elementary shellings (2^, deforming the boundary triangulation dA 
with induced action in the bulk. 

a. Bistellar moves. There are four bistellar moves in three dimensions: the (1,4), the (2,3) and their 
inverses. In the first, one creates four tetrahedra out of one by placing a point p in the interior of the original 
tetrahedron whose vertices are labeled pi, i = 1, ...,4, and by adding the four edges (p,pi), the six triangles 
(p,Pi,Pj)i^j, and the four tetrahedra (p,Pi,Pj,Pk)i^j^k- The (2,3) move consists in the splitting of two 
tetrahedra into three : one replaces two tetrahedra {u,p\,p2,P3) and {d,pi,p2,Pz) (u and d respectively 
refer to 'up' and 'down') glued along the (j>i,P2,P3) triangle with the three tetrahedra (u,d,Pi,Pj)i^j. The 
dual moves, that is the associated moves in the dual triangulation, follow immediately. See FIG. 2. 



(4,1) (3,2) 




FIG. 2: The (4, 1) and (3, 2) bistellar moves. 



b. Elementary shellings. Since the manifold (S^, T^) has non-empty boundary, extra topological trans- 
formations have to be taken into account to prove discretization independence. These operations, called 
elementary shellings, involve the cancellation of one 3-simplex at a time in a given triangulation (A, dA). In 
order to be deleted, the tetrahedron must have some of its two-dimensional faces lying in the boundary dA. 
The idea is to remove three-simplices admitting boundary components such that the boundary triangulation 
admits, as new triangles after the move, the faces along which the given tetrahedron was glued to the bulk 
simplices. Moreover, for each elementary shelling there exists an inverse move which corresponds to the 
attachment of a new three-simplex to a suitable component in dA. These moves correspond to bistellar 
moves on the boundary dA and there are accordingly three distinct moves for a three-manifold with bound- 
ary, the (3, 1), its inverse and the (2, 2) shellings, where the numbers (p, q) here correspond to the number 
of two-simplices of a given tetrahedron lying on the boundary triangulation. In the first, one considers a 
tetrahedron admitting three faces lying in dA and erases it such that the remaining boundary component 
is the unique triangle which did not belong to the boundary before the move. The inverse move follows 
immediately. The (2, 2) shelling consists in removing a three-simplex intersecting the boundary along two 
of its triangles such that, after the move, dA contains the two remaining faces of the given tetrahedron. 
These shellings and the associated boundary bistellars are depicted in FIG. 3. 
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(3,1) (2,2) 




FIG. 3: The (3, 1) and (2, 2) shellings, their dual moves and the associated boundary bistellars. 

The subset of bistellar moves and shellings which map T-adapted triangulations into T-adapted triangu- 
lations will be called adapted Pachner moves, and, considering a local simplicial structure 7^ = U„ =1 A3, 
k = 1, ...,4, an adapted (p, q) Pachner move T p i— > T q will be noted 3 3 (p, ( j), or more generally 3 . Any two 
PL-homeomorphic, T-adapted triangulations (A, dA) of the PL-pair T v ) are related by a finite sequence 
of such moves. 

Here it is important to take into account the PL-embedding of the string spin network graph T in the 
(dual) triangulation (A, dA). We will call = T n T£ the restrictions of the graph T to the local simplex 
configurations 7^ appearing in the moves. If the graph T^. is not the null graph, we will consider that it is 
open and does not contain any loop. If this was not the case, the set of adapted moves would reduce to the 
identity move, under which the transition amplitudes are obviously invariant. Hence, T^, if at all, can only 
be either an edge, or more generally a collection of edges, either a (three-valent) vertex. The associated 
string spin network functional ^r k will be represented by a group function <pk, which is the constant map 
0fe = 1 if Tfc is the null graph. We will make sure to check that, under a (p, q) Pachner move, 4>k transforms 
as 4> p h- > 4> q . 

2. Maximal tree moves 

It is also necessary to define topological moves for the trees [131 , [III ■ Any two homologous trees T\ and 
T2 are related by a finite sequence of the following elementary tree moves 7 : T\ \— ► T2. 

Definition 1 (Tree move) Considering a vertex Ao belonging to a tree T, choose a pair of edges A\, Aj 
in A touching the vertex A such that Ai is in T , A' x is not in T and such that A[ combined to the other 
edges of T does not form a loop. The move 7 consists in erasing the edge A± from T and replacing it by 




There is another operation on trees that we need to define. When acting on the simplicial complex (A, dA) 
with a bistellar move or a shelling, one can possibly map (A, cA) into a simplicial complex (A', dA') with 
a different number of vertices. Hence, a maximal tree T of A is not necessarily a maximal tree of A'; the 
Pachner moves have a residual action on the trees. This leads us to define the notion of maximal tree 
extension (or reduction) accompanying Pachner moves modifying the number of vertices of the associated 
simplicial complex. 
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Definition 2 (Tree extension or reduction) An extended (or reduced) tree T, associated to a Pachner 
move I 3 : A i— > A' modifying locally the number of vertices of a simplicial complex A, is a maximal tree of 
A' obtained from a maximal tree T of A by adding (or removing) the appropriate number of edges to T as 
a mean to transform T into a maximal tree Tj> of A' . 

Obviously, there is an ambiguity in the operation of tree extension or reduction. But, because of the fact 
that the regularized physical inner product will turn out to be independent of a choice of maximal tree, 
there will be no trace of this ambiguity in the computations of the transition amplitudes. 

3. Curve moves 

Finally, we define the PL analogue of the Rcidcmcister moves, which where in fact a crucial ingredient in 
the proof of Reidemeister's theorem. Any two ambient isotopic PL embeddings 71 and 72 of a curve 7 in 
the dual complex (A*, dA*) are related by a finite sequence of the following elementary topological moves 
g : 7j (-> 72. 

Definition 3 (Curve move) Consider a PL path 7 lying along the p boundary edges e\, e%, e p of a two- 
cell f of the dual pair (A*,9A*), where f has no other edges nor vertices traversed by the curve 7. Erase 
the path 7 along the edges e±, e%, e p and add a new curve along the complementary df \ {ei, €2, e p } of 
the erased segment in f. 



df\e 



We have now defined all of the elementary regulator moves. To summarize, an elementary reg- 
ulator move 51 : TZ(n,t) 1— > 1Z', e ,-> is a finite combination of all the above moves : !K[7^( 7)i£ )] = 
(T v , y[(A, dA)], 7[Tj>], G[a], C[C]). Proving the invariance of the regularized physical inner product (|45|) 
under all of these elementary regulator moves is equivalent to showing the independence of the regulating 
structure TZi^a. Note however that we have not included contractions or dilatations of the string tube T v 
radius r\ in the regulator moves. This is because the invariance under shellings implies the invariance under 
increasing or decreasing of 77. Indeed, the bistellars and shellings are the simplicial analogues of the action 
of the homcomorphisms Homeo[(E j; , T^)]. In particular, the topological group Homeo [(£,,, T,,)] contains 
transformations deforming continuously the boundary T v , like for instance maps decreasing or increasing 
the (non-contractiblc) radius 77 > of the boundary torus T^. Hence, showing the invariance under ele- 
mentary shellings is sufficient to prove the independence on the string thinckening radius, and the moves 
defined above are sufficient to conclude on the regulator independence of the definition of the regularized 
physical inner product. 

To push the result further and conclude on the topological invariance of the transition amplitudes, we 
need extra ingredients that we define in the following section. 



B. String spin network graph moves 



We now introduce the following elementary moves respectively acting on the edges and vertices of the 
open graph T. All ambient isotopic PL embeddings of the one complex T are related by a finite sequence of 
the following elementary moves noted S. 

Definition 4 (Edge move) An edge move is a curve move applied to an edge er of the graph T. 

Note that these moves apply also to the open edges of the graph V. However, there exists other moves 
which displace the endpoints. 

Definition 5 (Endpoint move) Considering an open string spin network edge er ending on the point 
Xk E X supported by a dual vertex v of dA* , which is such that its neighbouring vertex v' not touched by 
er belongs to dA* , an endpoint move consist in adding a section to er connecting v to v' . 
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We also need similar moves for the vertices. 

Definition 6 (Vertex translation) Let vy denote a three-valent spin network vertex sitting on the vertex 
v of the dual complex (A*, <9A*). Choose one edge er among the three edges emerging from vr and call v' the 
dual vertex adjacent to v which is traversed by er- Call e,e' C (A*,dA*) the dual edges locally supporting 
er, i.e, such that de = {v,v'} and v' = er H (e fl e'). 

The move consists in translating the vertex Vr along e from v to v' . This is achieved by choosing one dual 
face sharing the dual edge e and not containing the dual edge e' , and acting upon it with the edge move. 




Note that the use of rectangular faces in the above picture is only for the clarity of the picture, the move 
is defined for faces of arbitrary shape. 

It is important to remark that the above moves respect the topological structure of the embedding 
because no discontinuous transformations arc allowed and because the number and nature of the crossings 
are preserved since the faces used to define the moves are required to have empty intersections with the 
string or graph a part from the specified ones. 

The combination of the adapted Pachner moves and the spin network moves are the simplicial analogues 
of the action of the homeomorphisms Homeo[(£ r) , T v )] on the triple ((E^, T^), F). 

C. Invariance theorem 

We can now prove the following theorem. 

Theorem 1 (Invariance theorem) Let ^r denote a string spin network element of a given basis ofTLkin 
defined with respect to the one-complex T. Choose a regulator TZ(n,e) = (T ri ,(A £ ,dA £ ),T,a,C) consisting 
in a thickening T n of the string, a cellular decomposition (A e ,<9A c ) of the manifold (E^X^) adapted to the 
graph T, a maximal tree T of A, a closed path a in OA*, and a collection C of open paths (3 in dA* . Let 
H : ~R-(ri.t) l— * ^'(n e') ( res P- S : -T > V ) denote an elementary regulator move (resp. a string spin network 
move). The evaluated linear form (|46[) is invariant under the action ofH and 9-' 

P[K M -^ T ] = P[%[TZ M };^r] (50) 
= P[ft ( „ i6) ;* g( D]. 

Proof. We proceed by separatly showing the invariance under each elementary regulator moves, before 
proving the invariance under graph moves. 

• Invariance under maximal tree moves. 

Here, we simply apply to the proof of invariance under maximal tree moves writen in (26j . Firstly, we 
need to endow the tree T of the left hand side of the move with a partial order. To this aim, we pick 
a distinguished vertex r of T, chosen to be the other vertex of the edge A' : . The rooted tree (T,r) 
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thus acquires a partial order ^ : a vertex A' of (T, r) is under a vertex A , A < A , if it lies on the 
unique path connecting r to Ao- We can now define the tree Ta to be the subgraph of T connecting 
all the vertices above Ao : Ao is the root of Ta . The second ingredient that we need is the notion of 
Bianchi identity (|4"4")l applied to trees. Indeed, if T is a tree of the (regular) simplicial complex A, its 
tubular neighborhood has the topology of a 3-ball and its boundary has the topology of a 2-sphcrc. 
This surface S can be built as the union of the faces / dual to the edges Ai in A touching the vertices 
of T without belonging to T. Hence, applying the Bianchi identity to the tree Tj\ yields 

9f&=(ll9f)gfa (Usf). (51) 
\/eSi / \fes 2 ) 

Here, /o and f' Q are the faces dual to the segments Ao and A (note that Ao does not belong to Ta ). 
The sets <Si, S2 are the set of faces dual to the segments Ai touching the vertices of Ta without 
belonging to Ta , and which are not f nor /q. The presence of two different sets Si and S2 is simply 
to take into account the arbitrary positioning of the group element g f amoung the product over all 
faces. As usual, the group elements are defined up to orientation and conjugation. Next, we apply a 
delta function to both sides of the above equation and multiply the result as follows : 



n s (9f) = n s 

f g T f g T 

f # fo . /q / # /o > Jo 

n w = n 6 (9f)- 

f$T f<£T' 



n 9f 9f n 9f 

,/e-Si / V/GS2 



(52) 



In the second step, we have simply used the delta functions with which the expression has been 
multiplied to set the group elements associated to the sets Si and 1S2 to the identity (the faces of 
Si and S2 are dual to segments not belonging to T). One can then check that the various steps of 
the proof remain valid if the boundaries of the dual faces carry string spin networks. This shows the 
invariance of the regularized inner product (|46p under maximal tree move □ 

Invariance under adapted Pachner moves. 

To prove the invariance under Pachner moves, we introduce a simplifying lemma [23[, [24[. 

Lemma 1 (Gauge fixing identity) To each vertex of the dual triangulation A* are associated four 
group elements {g }a=i,...,4, six unitary irreducible representations {j a b}a<b=i,..A of G and a string 

spin network function (f>i({g a } a= i 4). If 4>i is the constant map <pi = 1> or depends on its group 

arguments only through monomial combinations {g a gb}a^b of degree two, then the following identity 
holds : 

4 f jab 4 [ jab 

Yl / dg a 3 ir (g a gb) &i({gd}d) = Yi / d 9aS{g c ) J ^ (g a gb) 4>i({g<i}d), (53) 
for c = 1, 2, 3 or 4. 



Proof of Lemma 1 . The above equality is trivially proven by using the invariance of the Haar measure 
and performing the change of variables j co = g c gb, for c < b (resp. 7& c = gbg c , for c > b) in the left 
hand side. This translation is always possible since the group function <fi is either the constant map 
or depends on the group elements only through monomials of degree two □ 

Let us comment here on the validity of the hypothesis made on the spin network function <f>i associated 
to the graph Ti = r D Ti, with Ti = A3 in the above Lemma. In fact <pi depends necessarily on 
combinations of the form g a gb locally if the graph Ti is not the null graph. Indeed, Ti can either be 
a collection of edges, in which case this requirement simply states that the edges are open, either a 
vertex, where one can always use the invariance of the associated intertwining operator to satisfy the 
desired assumption. Hence, this requirement is always locally satisfied. 

We can now show the invariance under bistellar moves and shcllings. 
— Bistcllars moves: 
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* The (4,1) move. Consider the four simplices configuration T4 in (A,<9A) (FIG. 2). Since 
the amplitudes do not depend on the maximal tree T of A, we are free to chose it. The 
simplest choice consists in considering a maximal tree T whose intersection T4 with the 
simplex configuration T4 reduces to the four external vertices and to a single one-simplex 
touching the central vertex. We work in the dual picture and label the four external dual 
edges from one to four. The face dual to the internal tree segment is chosen to be the face 
142. We note g a and h a b, a,b = 1, 4, the group elements 11 associated to the external and 
internal dual edges respectively, while the representations assigned to the dual faces are noted 
j a t,. The general PL string spin network state restricted to the configuration T4* is noted 
04({3a}a, {h a b} a<b) ■ Obviously, <f>i is not a function of all of its ten arguments, otherwise it 
would contain a loop, but can generally depend on any one of these ten group elements, as 
suggested by the notation. The regularized physical inner product (|4"6"|) restricted to these 
four simplices yields 

/ dgidg2dgy,dgi dhi2dh\ 3 dhi4dh2 3 dh2idh 3 i (54) 
Jg 10 

jl2 jl3 jl4 J23 324 ^34 

t (.91^1232) 7T (gih 13 g 3 ) 7T {g 1 h 14 g 4 ) n (g 2 h 23 g 3 ) tt (g 2 h 2 4g4) ^ (33^3434) 
S(hi2h 2 3h 3 i) 5{hi 3 h 34: h4i) S(h2 3 h 34 h 2 4) (j>i{{g a }a, {Kb}a<b), 

where we have omitted the sum over representations weighted by the associated dimensions. 
We start by implementing Lemma 1 at vertices 2, 3 and 4 to eliminate the three variables 
hu, b 7^ 1. We obtain 

/ dg 1 dg2dg 3 dg 4 dh23dh24dh 3 4 (55) 
Jg 7 

3 k (3132) J 7r 3 (.9133) J tt 4 (.9134) J ™ (32^2333) (32^2434) 3 k (33/13434) 
6(h 2 3) S(h 3 4) 5(h 23 h 3i h2i) 4>{{g a }a, {^23, ^34, h^}) 

= i dg 1 dg 2 dg 3 dg i J n {gig 2 ) 3 n (gigs) ^ (3134) 
Jg 4 

J K (3233) ^ (3234) ^ (3334) 0l({3a}a), 

where we have integrated over the delta functions to eliminate the interior variables in 
the second step. The right hand side of the above equality corresponds to the one sim- 
plex configuration of the (4, 1) move with the associated maximal tree reduction (the ob- 
vious removal of the internal tree segment; T\ is given by the four vertices of the resulting 
tetrahedron). In other words, we have just proved the invariance under the transformation 
%!) : (T 4 , T 4 ) 1— » (7i, T[), with T k =TnT k and T = Ty (41) . 

* The (3, 2) move. Here, we consider the three simplices configuration T 3 (see FIG. 2) and 
chose a tree T intersecting T 3 only on its five vertices. Concentrating on the dual graph, we 
label the three vertices from one to three and respectively note 3", h a b, and J^f , a, b = 1, 3, 
a, (3 = 1,2, the external and internal group elements, and the representation labels. The 
associated string spin network state is called 4> 3 . The transition amplitude, restricted to these 
three simplices, yields (omitting the sum over representations and associated dimensions) 



/ dg\dg\dg\dg\dg\dg\ dh 12 dh 13 dh 23 (56) 
Jg 9 

,11 -'11 ,11 .'22 -22 -22 

J12 1 1 , 1 v J13 / 1 , 1 \ J23 / 1 , 1 \ J12 / 2 t 2\ Jl3 / 1 7 2\ J23 / 2 7 2\ 

k (3i«i23 2 ) n (3i«i33 3 ) n (32^2333) n (3i«i23 2 ) n (32^1333) n (3 2 «233 3 ) 

^ (dWl) ^ (9292) ^ (3333) S(h 12 h 23 h 13 ) (j) 3 ({ga}a,{ h ab}a<b)- 

Using the gauge fixing identity at the vertices 2 and 3 to eliminate the variables h\ ai a =/= I, 



The notations takes into account the appropriate orientations. 
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and solving for the delta function leads to 

/ dg\dg\dg\dg\dg\dg\ (57) 
Jg 6 

-•11 ,-22 -22 -22 

J12 / 1 l x J13 /I 1 x J23 /I 1\ J12 / 9 9\ ^13 / 2 2\ J 23 / 9 9\ 

^ (5i5 2 ) t (.9i5 3 ) t (ff 2 3 3 ) t (3i5 2 ) n (3 2 3 3 ) ^ (ftjSs) 

^■12 jl2 jl2 

* (dial) % {g\g\) ¥ {glgl) <j>({g a } a ) 
dg\dg\dg\dg\dg\dg\ dh 



0l2 , 1 h i 1 In ^23 / 1 1\ J 12 / 2 2\ J 13 / 2 2\ •' 23 / 2 2\ 

Jll J22 J 12 

^ (fl£%i) % {g\hgl) tt {g\hgl) M{9a}a,h), 

where we have used the inverse gauge fixing identity in the last step. This expression corre- 
sponds to the two simplices configuration T 2 of the (3, 2) move. 

— Shcllings : 

* The (3, 1) move. Remarkably, writting the amplitudes associated to the left and right hand 
sides of the (3, 1) shelling leads to the same expression than the (4, 1) bistellar, even if the 
geometrical interpretation is obviously different. This is due to the fact that we are imposing 
the flatness constraint F = also on the faces of the boundary 12 simplicial complex <9A and 
integrating also on the boundary edges. The only difference is in the presence of possible 
open string spin network edges reflected in the group function = (f>({g a } a , {A a }a), without 
any incidence an any steps of the proof given for the (4, 1) bistellar. Accordingly, the proof 
of invariance under the (3, 1) shelling is the one sketched above. 

* The (2, 2) move. The same remark applies here, the amplitudes are exactly identical to the 
ones of the (3, 2) bistellar. 

Accordingly, we have proven the invariance under adapted Pachner moves □ 

• Invariance under curve moves. We here show that the regularized physical inner product is invariant 
under curve moves. The proof uses the flatness constraint F = 0. Consider a particular dual face / of 
(A*,<9A*) containing n boundary edges positively oriented from vertex 1 to vertex n. Suppose that 
there are p < n dual edges e±, e p supporting a curve positively oriented w.r.t the face /, to which 
is associated a spin j representation. We want to prove that the associated amplitude is equal to the 
amplitude corresponding to the curve lying along the n — p edges of df \ {e%, ...,e p } after the edge 
move. We start from the initial configuration 

c n 

/ II dg a tt (g\...g p ) S(gi...g p g p+ -i...g n )5(g 1 Gi)S(g 1 H 1 ) ... 8(g n G n ) S(g n H n ), (58) 

where the capital letter G a , H a , a = 1, n, represent the sequences of group elements associated to 
the two others faces sharing the edge a. We then simply integrate over the group element g\ to obtain 

- Tl 

/ I dg a i\ (g- 1 ..^) Sig-K-.g^G,) S(g^ l 1 ...g^ 1 Hi) ... 5(g n G n ) 5(g n H n ) (59) 
= ~[dg a tr (g^.-.g'^) <%i-. ) 5( 9l Hi) ... 5{g n G n ) 5{g n H n ). 

•^ G " o=l 

Note the reversal of orientations intrinsic to the move. This closes the proof of invariance under curve 
move □ 

We finish the proof of theorem 1 by showing the second part, namely the invariance of the transition 
amplitudes under string spin network graph moves. 



this sense, the boundary amplitudes are very different from a 2 + 1 quantum gravity model defined on an open manifold. 
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• Invariance under edge moves. 
The proof is the one given for the curve move □ 

Invariance under endpoint moves. 

Here, we use the momentum conservation Dp = 0. Considering a particular dual face / containing 
n boundary edges, with p < n dual edges ei,...,e p supporting an open string spin network edge 
(positively oriented w.r.t /) ending on the boundary of the edge p, to which is associated a spin j 
representation. We call the string field evaluated at the target of the kth edge. We choose the 
holonomy starting point x to be on the endpoint of the p edge (we prove below that nothing depends 
on this choice) and, since nothing depends on the paths (3 by virtue of the invariance under curve 
moves, we choose a path (3 of C along the edge p + 1. The relevant amplitude is given by 

„ Tl 

/ Y\_dg a d\ p d\ p+ i n (g 1 ...g p \ p )S(g p+1 Xp + i\p 1 ), (60) 
U ' G " a=l 

where the notations are the same than above. It is immediate to rewrite the above quantity as 

J dg a d\ p d\ p+1 tt {g 1 ...g p g p+1 X p+ i)5{g p+ i\ p+ i \~ : ), (61) 

J " a=l 

which concludes the proof of endpoint move invariance □ 

• Invariance under vertex translations. 

Here, we consider three dual face fi, i = 1,2,3, of (A*,i9A*) each containing ni boundary edges 
positively oriented from vertex 1 to vertex n^. The three faces meet on the common edge e which is 
such that 1 = t(e), i.e., e = e\., forall i. Suppose that there arc p\ < n\ dual edges e\, •..,e^ i (rcsp. 
Pi < n 2 dual edges e\ , e^ 2 ) of the face /i (rcsp. f 2 ) supporting a string spin network edge ep (rcsp. 
ep) colored by a spin ji (resp. j?) representation and oriented negatively w.r.t the orientation of fx 
(resp. fa)- Suppose also that the face contains 713 — p% dual edges e^ 3+l , e along which lies a 
positively oriented (w.r.t. to the orientation of the face) string spin network edge ep colored by a spin 
j'3 representation. Consider that the three edges meet on the vertex vr supported by the vertex 1 
of (A*,9A*). Noting g the group element associated to the common edge e, one can write the spin 
network function associated to the three valent vertex vr lying on 1 and use the invariance property 
of the associated intertwining operator 1 to 'slide' the vertex along the edge e : 

n \9p 3 +i~-9n 3 -i9) hifoh 

It is then possible to use the flatness constraint F = on either of the faces /1 or f% to implement an 
edge move on ef or ep thus completing the vertex move □ 

By virtue of all the above derivations, we have now fully proven theorem 1. 

To be perfectly complete, we need to verify that the amplitudes are also independent under orientation 
and holonomy base point change. This leads to the following proposition. 

Proposition 1 The regularized physical inner product (|46[) is independent of the choice of orientations of 
the dual edges and faces o/(A*,9A*), and does not depend on the choice of holonomy base points. 

Proof of Proposition 1. It is immediate to see that the amplitudes do not depend on the orientations 
of the dual faces and dual edges of (A*<9A*), nor on the holonomy starting points on the boundaries of 
the dual faces [26[. Indeed, a dual face and dual edge orientation change correspond respectively to a 
change g e 1— > g~ 1 and gf 1— > gj 1 which are respectively compensated by the invariance of the Haar measure, 

dg e = dg~ x , and of the delta function: S(gj) = S(gJ 1 ). A change in the holonomy base point associated 
to a dual face / will have as a consequence the conjugation of the group element gf by some element h in 
G. Since the delta function is central, 8{gf) = Sfigjh' 1 ), the regularized physical inner product (|4"6")) will 
remain unchanged under a such transformation. 
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Concerning the base point x used to define the holonomics along the loop a and the paths f3 G C in (|46[) . 
the situation is similar. Let x be noted x\ and suppose that we change the point X\ to another point X2 
in X neighbouring x\. Since we have showen the invariance under bistellars and shcllings, we are free to 
choose the simplest discretization 13 of the manifold T r; ). We choose it such that the cylindrical section 
of between x\ and xi is discrctized by a single dual face with two opposite sides glued along a dual edge 
e linking x\ to xi. By virtue of the curve move invariance, we are also free to choose the path (3 to be along 
e. The amplitude based on X2 as a starting point for the paths a and /3, restricted to this section of T^, 
yields 

f 2 

/ r II d 9ad\adg,3 5{g2\2uX 2 1 )5{gp\iX 2 1 )5{g 2 g,39i 1 9j i 1 ) f{{ga\a,{^a\a,gp), (63) 

J* 35 a=l 

where g a is the holonomy around the disk bounding the tube section at the point x a and the function / 
describes the string spin network function together with the other delta functions containing the group 
elements g a and g@. It is immediate to rewrite the above expression as 

/ r n d 9adXadg(3 5(g 1 \ 1 u\ 1 1 )6(g f3 1 \ 2 \ 1 1 )8(g 1 1 g l3 1 g2gp) f({g a }a,{^a}a,gp), (64) 

a=l 

which is the amplitude based on x\ as a starting point for the paths a and [3 □ 

There are two major consequences due to the above theorem and proposition. First, there is no continuum 
limit to be taken in (|45|) . Since the transition amplitudes are invariant under elementary regulator moves, 
the regularized physical inner product ([46]) is independent of the regulator and the expression (|46|) is 
consequently exact, there is no need to take the limits 14 e, T) — > 0. In particular, we have shown that 
the amplitudes are invariant under any finite sequence of bistellar moves and shellings which implies, by 
Pachner's theorem, that the physical inner product is well defined and invariant on the equivalence classes of 
PL-manifolds 15 (A, <9A) up to PL-homeomorphisms. Accordingly, the transition amplitudes are invariant 
under triangulation change and thus under refinement. This leads to the second substantial consequence of 
theorem 1. The crucial point is that the equivalence classes of PL-manifolds up to PL-homeomorphism arc 
in one-to-one correspondence with those of topological manifolds up to homcomorphism. Sec for instance 
25| for details. Hence, showing the invariance of the regularized physical inner product under triangulation 
change is equivalent to showing homcomorphism invariance: the discrctized expression (|46j) is in fact a 
topological invariant of the manifold (£,,, T v ). In particular, the amplitudes are invariant on the equivalence 
classes of boundary torii T v up to homeomorphisms. It follows that they do not depend on the embedding 
of the string 5? . 

Combining these results with the second part of theorem 1 stating that the regularized physical inner 
product is invariant under string spin nctwotk graph moves, we obtain the following corollary. 

Corollary 1 The physical inner product (|46[) is a topological invariant of the triple ((S^,T r; ),r) : 

P[ft fee) ;*r]=P[[(£^)];* [r] ], (65) 



where [(£,,, X^)] and [T] denote the equivalence classes of topological open manifolds and one-complexes up 
to homeomorphisms and ambient isotopy respectively. 

This corollary concludes our study of the topological invariance of the theory of extended matter coupled 
to BF theory studied in this paper. 



13 Anticipating on the next paragraph, we are using the fact that the invariance under Pachner moves implies the topological 
invariance of the amplitudes. Accordingly, we can use a cellular decomposition which is not necessarily a triangulation. 

14 More precisely, we have shown that H46H is invariant when going from to IZf^i^i) for all (r)' ,e') ^ (r/, e) which implies 
regulator independence. 

15 More precisely, a triangulation A is not a PL-manifold. It is a combinatorial manifold which is PL-isomorphic to a PL- 
manifold. 
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V. CONCLUSION 



In the first part of this paper we have studied the geometrical interpretation of the solutions of the BF 
theory with string-like conical defects. We showed the link between solutions of our theory and solutions 
of general relativity of the cosmic string type. We provided a complete geometrical interpretation of the 
classical string solutions and explained (by analyzing the multiple strings solution) how the presence of 
strings at different locations induces torsion. In turn torsion can in principle be used to define localization 
in the theory. 

We have achieved the full background independent quantization of the theory introduced in [l[ . We showed 
that the implementation of the dynamical constraints at the quantum level require the introduction of 
regulators. These regulators are defined by (suitable but otherwise arbitrary) space discretization. Physical 
amplitudes are independent of the ambiguities associated to the way this regulator is introduced and are 
hence well defined. There are other rcgularization ambiguities arising in the quantization process that have 
not been explicitly treated here. For an account of these as well as for a proof that these have no effect on 
physical amplitudes see [3fj| . 

The results of this work can be applied to the more general type of models introduced in [3l[, were it 
is shown that a variety of physically interesting 2-dimensional field theories can be coupled to the string 
world sheet in an consistent manner. An interesting example is the one where in addition to the degrees of 
freedom described here, the world sheet carries Yang-Mills excitations. 

There is an intriguing connection between this type of topological theories and certain field theories in 
the 2+1 gravity plus particles case. One would expect a similar connection to exist in this case. However, 
due to the higher dimensional character of the excitations in this model this relationship allows for the 
inclusion of more general structures: only spin and mass is allowed in 2+1 dimensions. The study of the 
case involving Yang-Mills world-sheet degrees of freedom is of special interest. This work provides the basis 
for the computation of amplitudes in the topological theory. A clear understanding of the properties of 
string transition amplitudes should shed light on the eventual relation with field theories with infinitely 
many degrees of freedom. 
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